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Abstract 

The photon emission intensity spectrum is calculated taking into ac- 
count influence of multiple scattering (the LPM effect) under conditions 
of recent CERN SPS experiment. It is shown that the theory quite satis- 
factory describes data. The integral characteristics: the radiation length 
and total probability of photon emission are analyzed under influence of 
the LPM effect, and the asymptotic expansions of these characteristics are 
derived. 
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1 Introduction 



When a charged particle is moving in a medium it scatters on atoms. With 
probability ~ a this scattering is accompanied by a radiation. At high energy 
the radiation process occurs over a rather long distance, known as the formation 
length l c (see e.g. [I]): 

1 + m z c<j 

where u is the energy of emitted photon, e{m) is the energy (the mass) of a 
particle, e' = e — oj, d c is the characteristic angle of photon emission, the system 
h — c — 1 is used. 

Landau and Pomeranchuk were the first who showed that if the formation 
length of bremsstrahlung becomes comparable to the distance over which the 
multiple scattering becomes important, the bremsstrahlung will be suppressed 
j2J- Migdal jSj developed the quantitative theory of this phenomenon. 

A new interest to the theory of the LPM effect is connected with a very 
successful series of experiments performed at SLAC jl] . In these experiments the 
cross section of the bremsstrahlung of soft photons with energy from 200 keV to 
500 MeV from electrons with energy 8 GeV and 25 GeV is measured with an 
accuracy of the order of a few percent. Both LPM and dielectric suppression are 
observed and investigated. The logarithmic binning of photon spectrum was used 
with 25 bins per decade. These experiments were the challenge for the theory since 
in all the mentioned papers calculations are performed to logarithmic accuracy 
which is not enough for description of the new experiment. The contribution of 
the Coulomb corrections (at least for heavy elements) is larger than experimental 
errors and these corrections should be taken into account. 

Recently new study of LPM effect at higher energies of electrons (e =149, 207 
and 287 GeV), where the effect has influence upon much wider part of spectrum 
comparing with e =25 GeV, was performed in the H2 beam line of the CERN 
SPS 0,0- The logarithmic binning of photon spectrum was used with 25 bins 
per decade as at SLAC. 

We developed the new approach to the theory of the Landau-Pomeranchuk- 
Migdal (LPM) effect [7]. In this paper the cross section of the bremsstrahlung 
process in the photon energies region where the influence of the LPM is very 
strong was calculated with a term oc 1/L , where L is characteristic logarithm 
of the problem, and with the Coulomb corrections taken into account. In the 
photon energy region, where the LPM effect is "turned off", the obtained cross 
section gives the exact Bethe-Maximon cross section (within power accuracy and 
with the Coulomb corrections). This important feature was absent in the previ- 
ous calculations. The LPM effect in a thin target were the interference effects 
and boundary photon emission is of special interest was analyzed in jS]. The 
probability of multiphoton emission is enhanced at high energy jU] . Correspond- 
ingly this effect is very important for comparison of theory prediction and data. 
The influence of LPM effect on integral characteristics of bremsstrahlung was 
considered in [TO]. The other approaches to the LPM effect theory see e.g. in 
|llj . |12j . the recent review is given in |13j . 

In Sec. 2 the theory predictions are compared with the recent CERN SPS data 
It is shown that the theory quite satisfactory describes data. In Sec. 3 
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the variation of the radiation length due to multiple scattering is discussed. In 
the region where the LPM effect is weak (e <C e e ) the asymptotic expressions for 
the radiation length and the photon emission probability are derived taking into 
account both decomposition over e/e e and over 1/L\, L\ is the characteristic 
logarithm of the problem. Appendices A and B contains details of derivation 
for radiation length and Appendix C for the photon emission probability. In 
Appendix D the structure of series over e/e e and over 1/ L\ is analyzed for e <C e e . 



2 Influence of the multiple scattering on the 
bremsstrahlung spectrum 

The spectral radiation intensity obtained in [7j, Eq.(2.39) (see also Sec. 3 in [TUj ) 
is valid for any energy and has the form 
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where z — vt, p — i/(2v), ip(x) is the logarithmic derivative of the gamma 
function. The functions fi(z) and f2(z) are defined by following expressions 
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here Li 2 (x) is the Euler dilogarithm. Use of the last representation of function 
G(z) simplifies the numerical calculation. The crucial parameter v is 
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where Z is the charge of the nucleus, n a is the number density of atoms in the 
medium, A c = 1/m = (h/mc) is the electron Compton wavelength. Here e e 
is a characteristic parameter of medium, starting from this energy the multiple 
scattering distorts the whole spectrum of bremsstrahlung including its hard part: 
for iridium e e =2.27 TeV, for gold e e =2.6 TeV, for tungsten e e =2.73 TeV, for 
tantalum e e =3.18 TeV, for lead e e =4.38 TeV. 

In the case e<e e the LPM effect manifests itself when 

ux(x c ) = 1, x c = — ^— ~ — . (2.5) 

In this case in the hard part of spectrum (1 > x ^> x c , v\ ~ x c /x <C 1) one has 
(see Appendices A and B) 
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here L^Jf is the Bethe-Maximon radiation length. Note that if neglect here the 
terms oc e 2 /e 2 we obtain the Bethe-Maximon intensity spectrum. 

In the recent CERN SPS experiment the iridium target with thickness Z=0.128 mm 
(I I ' L rad = 4.36 ± 0.10%, L rad is the radiation length) [3] and tantalum target [S] 
with thickness l/L rad = 4.45 ± 0.10% were used. Photons with energy 2 GeV 
< ijj < e were detected in a lead glass calorimeter. 

Obtained experimental data should be recalculated: 
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Because photon energies were histogrammed logarithmically, using 25 bins per 
decade of energy, one has for the coefficient k 

k h = = e — 1 = 0.096, k m = 2 = 2— = 0.092, 

^min ^rnax ^min C s -\- L 

j ^max ^min 1 -s n oo 1 ^max l n 10 , _s 

ki = = 1 — e = 0.88, s = In = = 0.092, (2.8} 

^max Ldmin 25 

depending on the normalization point within bin. Here we will use k m . 

Since only high energy photons (u > 2 GeV) were measured, no boundary 
effects were observed. So, only the pure LPM effect was studied. For iridium 
e e =2.27 TeV and from one has Eq. (J2.5|) that the characteristic photon energy 
uj c (e) for which the LPM effect is well manifests itself is u c (287 GeV)=32 GeV and 
u c (207 GeV)=19 GeV, while for tantalum e e =3.18 TeV and w c (287 GeV)=26 GeV. 
The results of calculations for iridium (the initial electron energy e=287 GeV 
and e=207 GeV) and for tantalum (the initial electron energy e=287 GeV) are 
shown in Figs. 1-3. The curve 1 is the Bethe-Maximon intensity spectrum (see 
e.g.Eq. (|2.6jl ). The curves 2, 3 are calculated using Eqs. (|2.1j) and ()2.2j) . The 
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curve 2 presents the main term (the function <&(i>)), the curve 3 presents the 
correction term (the function F(v)). It should be noted that the prediction of 
our theory (sum of previous terms, curve 4) in the hard end of spectrum coincide 
with the Bethe-Maximon curve within the accuracy better than 10~ 3 . It is seen 
that given values uj c (e) show the scale where LPM effect becomes essential. For 
used thickness of target the multi-photon effects are very essential. The reduc- 
tion factor / in this case it is convenient to calculate using the following general 
expression (Eq.(3.4) of (5]): 
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where the main term (with $(v)) of spectrum Eq. (|2.1)l is substituted. The results 
obtained are in a good agreement with Eqs.(2.11) and (2.26) of [9J. The reduction 
factor / in iridium for three used energies is presented in Fig.4. It follows from 
Fig. 4 that if the electron energy decreases the reduction factor / diminishes as 
well. The final prediction, with the reduction factor taken into account for used 
target thicknesses, present the curves T. The data are recalculated according 
with Eqs. fl2.7j) . (j2.8|) using the coefficient k m . In Figs. 1-3 one can see that for 
energy 287 GeV there is quite satisfactory agreement of theory with data for 
both iridium and tantalum, for energy 207 GeV in iridium the agreement is 
somewhat less satisfactory. 



3 The LPM radiation length and the average 
number of emitted photons at e <C e e 

The local radiation length L ra d is defined by equation 

% = - J ( £ ) = -T- ^ 

at Wad 

In the absence of LPM effect the value L rad doesn't dependent on energy and 
defined by Eq. ()2.6|) . In this case, after averaging Eq. ()3.1jl with the electron 
distribution function over energy f(s,t) one obtains the closed equation for the 
mean energy loss which has the solution (e(t)) = eoexp(—t/L rac i). 

For the first photon emission according with Eq. (j3.1j) L rad = e/I(e). The 
dependence on the energy of the function (I(e)/e)L^ = L^J / L raa - for iridium 
(curve 1) and lead (curve 2) is given in Fig.5. The 1(e) is the integrated over x 
Eq. fl2.1|) . The relative value of correction (the term with F(v)) to the main term 
(with Q(v)) depends on energy. It attains the maximum ~ 5% near e ~ e e for 
both media. So, the value of L ra d for the energy e = 3 TeV increases 1.33 times 
in Ir and in 1.16 times in Pb, while for the energy e — 10 TeV the corresponding 
figures are 1.58 and 1.38. 

In the region where the LPM effect is weak (see Appendices A and B) 
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The first term of expansion over e/e e (not over 1/Li) was found in [TU] . 

The accuracy of Eq. (j3.2|) is defined by a several essentially different factors. 
The first is the relative value of discarded terms 0(e 3 /el) in Eqs. (jA.13|) . IB. 15]) . 
The second is the substitution of L c by L\ (g c = 1) in Eqs. (j2.1|) - (|2.4|l for the 
whole spectrum, the relative accuracy of this substitution is 0(e/ (e e L\)). The 
third is due to fact that in the initial formula Eq. (|2.1j) the terms oc 1/L% were 
rejected. However, as is shown in in Appendix D, in the region e <C e e the 
corrections oc 1/Lf are beginning with the terms contained e/e e . So, in this 
region, the relative accuracy in the formula Eq. (|2.1J) is also 0(e/(e e Lf)). 

One can see from Eq. ()3.2|l . that in a decomposition over e/e e one have to 
take into account the terms oc 1/Li because of relatively large numerical values 
of the coefficients of decomposition. Moreover, for heavy elements L\ ~ (6 + 7) 
taking into account of term oc 1/Li results in substantial change of numerical 
value of coefficients at the given degree of e/e e . Particularly this is important 
for term containing \n(e e /e). As a result, there is essential compensation inside 
the coefficients of decomposition. This permit to use the decomposition Eq. ()3.2|) 
at relatively large energies. At e/e e < 1/10 an error in the correction terms 
doesn't exceed 10%. 

As illustration we present the radiation length in iridium (Li = 6.9225, e e = 
2.27 TeV, L^ f (Ir)=2.91 mm): 



! ' ' 1 -0.464- + 0.045^ In - + 1.492^ J . (3.3) 



Zw(Ir) L™(Ir) \ e e ' e 

It should be noted that after substitution of Eq. ()3.3j) into Eq. (|3.1j) and averaging 
with the distribution function f(e,t) the higher moments of e appear. Because 
of this, equation ()3.1j) ceases to be closed. We will discuss this item elsewhere. 

The mean value of photons emitted by electron is also of evident interest. Par- 
ticularly, at uj <C uj c ~ e 2 /e e this value defines mainly the exponent in Eq. (J2.9|) . 
At energy where the LPM effect is rather weak (e e e ) the radiation probability 
per unit time integrated over whole photon spectrum is (see Appendix C) 
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It is seen that allowing for terms oc 1/Li changes essentially non-logarithmic 
terms in Eq. (J3.4|) . 
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A Appendix 

We consider here the asymptotic behavior of the total intensity of radiation J in 
the region e <C e e where the LPM effect is weak and one can put in Eqs. (j2.1j) - (|2.4)) 
v 2 = ivf, g c — 1, L c — L\. This substitution will manifest itself only starting from 

terms of the order — Iq. Then in expressions for $(z^) and F{v) in Eq. (j2.2|) we 

integrate by parts the terms oc r± and make the substitution of variable t — > —it. 
As a result we obtain 

: f Jm + ^J v y (A.l) 
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here the function G(z) is defined in Eq. (l2.3j) and 
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In Eq. ()A.2|) one can expand the function <px(z) up to terms cx z 4 . We find 
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In the integral with function tp 2 (z) in Eq. (|A.2|) we perform twice integration by 
parts over t. We find 
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To calculate S we choose the value y in such way that y <C 1 and \/3\y 3> 1 
and divide the integral over x into two. The first one is 
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Summing 5i, 5 2 i and 5 22 we obtain for S 
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Substituting Eq. (jA.ll|) into J m2 , adding with J m i Eq. ()A.5|) and taking into ac- 
count that 

^(0) = ^ 4) (o) = i y. 2 (o) = 4, ^ 2) (o) = A, ^ 2 4) (o) = H> 



15 



63 



oo 

J \nz(p i 2 ) {z)dz = -0.286, 



(A.12) 



we obtain 

, 4vr e QAe 2 (. e e 5 7 63 

= l_iLl + 5if!f.„£i- 2.040) + of 
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B Appendix 

Here we calculate the function J v in Eq. (|A.l|) defined in Eq. (|A.3|) . In the last 
expression one can expand the function li(z,t) up to terms cx z A . We find 
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where the function ip(x) is defined in Eq. (|2.2j) . The integral with the function 

h(z,t) is 
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The term with b 2 {z) can be calculated as the term J m2 in the Appendix A. 
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Now we turn over to the term with a 2 (z) in Eq. (jB.2 
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As in previous Appendix we divide the integral over x into two: < x < Xq and 
x <x <1. We have (see Eas. (lB~Tl) and (jA~H) ) 
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In the second integral we perform twice integration by parts over t and pass to 
the variable y. We find 
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We proceed as above (see Eq. ()A.8jl ) and expand (1 + y 2 ) 3 ~ 1 — 3y 2 , then 
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Integrating by parts over t two times we obtain 
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The term with the function g(z) can be calculated in the same manner as in 
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Taking into account that g(0) = — a^\o) and summing T\ Eq. ()B.5|) . T 2 \ Eq. ()B.7|) 
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Substituting the obtained expression for T Eq. flEU) and 4? Eq. ffiOf into J v 
adding to J v \ and taking into account that 
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C Appendix 

Here we consider the asymptotic behavior of the total probability of radiation 
W in the region £ C £ e to within terms ~ e/e e inclusively. According with 
Eqs. (lA.lJ) - (|A.4J) we present the total probability in form 
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The function entering in Eq. ([C.2|) are defined in Eqs. (|A.3|) - (|A.4|) . 

The term with (pi(z) in the integrand of W m to within mentioned accuracy is 
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In the integral with ^2(2) we divide the integral over x into two 

1 



4Re J — (1 - x) J dte- t ip 2 (z) ~ 4(lnx + l)<p 2 (0), (C.4) 
3^0 



and 



/5j/o , 00 



W^l = -8Re 



00 

+ 16Re y ydy y dte' 1 ^ 



yo 00 





,(2) 



(C.5) 



Integration of W„Ji coincides with calculation in Eq. (|A.10j) with substitution 



(4) 

¥>2 



^2 [ - ] ■ We have 



w {2) 

vv m21 



^ Re y ^ fdte-^A- 



00 

In j + 2C + lnyg) <^ 2 (0) - 2 y lnz^^tfc 



(C.6) 



Integrating by parts twice over t in integral W^ 2 Eq. ([TT5jl we find (see Eq.(fA~9f) 

/3yo , 00 



= 16-Im / ^ / dte-Vf - ) - 4vr^(0). 



(C.7) 





Summing En. (jU^|) . Eq.JCSI) and Ea.<|rT7|i and taking into 

account Eq. (|A.12|l ) we have 



W m = -4^ 2 (0) 



£ 4^(0) ^ 2 (0) 



+ 



4+ *?>«>) = i 

£ e O 
4 /, £ 



, £ e _ 5 /" , / 1 cosh z x 

In — + 2(7- - + 12 /In* — -— -3- 
£ 8 7 \ 2 J sinh z 



fi2 



87T £ 
15 £ fi 



- (In - + 1.959 ] + — - 

15 e f 



57T £ 



The first term in Eq. (|C.8|l ) is similar to Eq.(2.21) in 
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We will consider now the probability W v in Eq. (jC2|0 . The term with l\(z, t) 
gives no contribution into two first terms of expansion over — . 

Se 

1 oo 

W v2 = -4Re J {I- x)^ J dte- 1 [a 2 (z)(lnt + C) + b 2 (z)} = W$ + W$ '. (C.9) 
o o 

The calculation of the term W^ 2 (with 6 2 (z)) is analogous to the calculation 
above (with substitution (p 2 (z) — > b 2 (z)) 

oo 

W$ = -462(0) (iny + 2C-l) +8^ ln^ 1) (^)rfz + 47r^4 2) (0). (CIO) 
For calculation of term Wj; 2 (with 02(2)) it is convenient to use the variable 



wty = W { o\ + W ( ol 

' v v2 v2l ' v22 > 



00 , 00 



vv v2l 



WM = 8Re / ^ I 1 - 



^ /dte-«(a 2 (-J-a 2 {-) )(lnM C). 



V 

w 



v22 



dy 



y 



1 + 3/ 



2A2 



dte-\a 2 {z) -a 2 (0))(lnt + C). (C.ll) 



The integral over y in W^ 21 is the integral Frullani 



«2 



a 2 



,yJJ y 



-a 2 (0)lnt. 



(C.12) 



Using the last result in Eq. (|B.ll|l ) we have 



W41 = 3vr 2 a 2 (0). 



(C.13) 



In the expression for wj 22 we divide the integration interval into two parts: 



= 8Re 



WO 



dy 

y 



1 - 



i + y 



2^2 



/ 2 \ 

/^(4 2) (0)^ + ...)(lnt + C)^ = 0, 



yo 00 



W^l ] ~ 16Re y yrfy J e-\a 2 {z) - a 2 (0))(lnt + C)dt 







161m 



0yo , 00 
dy 

y 

y 



a 2 



(lnt + C) + d 



dt 



(C.14) 



The function d(z) is defined in Eq. flB.6j) . the further calculation can be done as 
in Eq.flHH) 



^ ~ 4vr-d(0) = Q7T-a ( 2 >(0). 



(C.15) 
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Summing W$ Eq. (|CT^ . Eq.JCH and W$\ Ea. ljTHSj) we obtain the 

following expression for W v Eq. fjC.lJ) 

W„ = -44 a (0) (in - + 2C - l)+8 f \nzb 2 1] (z)dz+—+2n-(3a { 2 2 \o)+2b { 2 2) (0)). 

\ S J J 3 Eg 

(C.16) 

Taking into account that 

00 

/ \nzb ( 2 \z)dz = 0.0431 


and using Eq. ()B.13|) we obtain (02(0) = 1/3) 

^ ln? + , 7655 _^£ + (|). (c , 7) 



D Appendix 

Here we discuss a structure of corrections to the total intensity of radiation at 
e <^ e e which is not included in the expansion over 1/L C contained in Eqs. (j2.1jl - 
(|2.4j) . We will show that these corrections, which are contained terms cx 1/L\ 
and higher degrees of 1/Li, are beginning from terms of the order e/e e . We will 
use the general expression Eq.(2.12) of 7] before decomposition over 1/L: 

I = (0|n (G- - G^) + r 2 p - G^) p| ) , (D.l) 



where 



el — x 2 ( T 1 £? 2 



G = p + 1 - iV(e), Go = + 1, 7(e) = — Q 2 [Li -ln^ -2C1 , 

(D.2) 

here the quantities r 1; r 2 are defined in Eq. ()2.2|) and e e and are defined in 
Eq. ([2.4)1 . In the region 1 > x ^> s/e e we expand the combination entering in 
Eq. (|D.1[) over degrees of V 

G- 1 - G^ = GqHVGq 1 + GqHVGoHVGq 1 + G^iVG^iVG^iVG^ 1 + ... 

(D.3) 



Substituting Eq. (jD.3|) into Eq. (jD.l|) one can verify that the terms of the order 
1/L 2 and higher degrees of 1/Li are appeared starting from the third term of 
decomposition Eq. (jD.3|) . Their relative contribution is of the order e 2 /e 2 . Taking 
into account that 

00 

G 1 = i [ e- u e- ip2t dt, < ei \e- lvH \Q 2 >= e *(ffi-ft) 2 , (D.4) 

J " Ait it 



we obtain 



1 °° 

< eil^o 1 !^ >= j- / ex P 




1 , , 2 



~ it + ~AV Bl ~ 82 



^1-K (\q 1 -q 2 \), (D.5) 
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where K (z) is the modified Bessel function (Mac-Donald's function). Substitut- 
ing this result into Eq. (jD.3|) and then into Eq. (jD.l|) we get 



dl am 2 x 



dx 7r(l — x) 



am 2 e fi + 2r 2 



47T S e 



nK 2 (g) + r 2 Kl(g)] V(g)gdg 
r 2 - r 1 



Li + 



9 



Obtained here result for agrees with Bethe-Maximon formula. 

In the region x ~ e/e e the potential V ~ 1 and decomposition Eq. (|D.3Jl is 
inapplicable. However in this region 



dl ~ am' 



oc doc 
1 — X 



~ am 



(D.7) 



and the relative contribution of this region ~ e/e e . 
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Figure captions 



• Fig.l The energy loss spectrum — in units — ^ — in iridium target with 

thickness Z=0.128 mm (//Lf a A / = 4.36%) for the initial electrons energy is 
£ = 286.6 GeV. The Coulomb corrections are included. 

— Curve 1 is the Bethe-Maximon spectrum, 

— curve 2 is the contribution of the main term describing LPM effect; 

— curve 3 is the correction term; 

— curve 4 is the sum of two previous contributions; 

— curve T is the final theory prediction with regard for the reduction 
factor (the multiphoton effects). 

Experimental data from [5], [U]. 

• Fig. 2 The same as in Fig.l but for the initial electrons energy e = 206.7 GeV . 

• Fig. 3 The same as in Fig.l but for tantalum target with thickness (Z/L^Jf = 
4.45%). 

• Fig. 4 The reduction factor for iridium target with thickness / = 0.128 mm 
(4.36 % L rad ) versus x = ui/e. The curves 1, 2, 3 are for energies 149 GeV, 
207 GeV and 287 GeV correspondingly. 

• Fig. 5 The relative energy losses of electron per unit time including the 
contribution of the correction term (1(e) /e) iff = L^/L ra d in iridium 
(curve 1) and in lead (curve 2) vs the initial energy of electron. 
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